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Cu2(l,4-diazacycloheptane)2Ci4 contains double chains of spin-i Cu 2+ ions. We report ac sus- 
ceptibility, specific heat, and inelastic neutron scattering measurements on this material. The mag- 
netic susceptibility, x{T), shows a rounded maximum at T — 8 K indicative of a low dimensional 
antiferromagnet with no zero field magnetic phase transition. We compare the x(T) data to exact di- 
agonalization results for various one dimensional spin Hamiltonians and find excellent agreement for 
a spin ladder with intra-rung coupling Ji = 1.143(3) meV and two mutually frustrating inter-rung 
interactions: J2 = 0.21(3) meV and J3 = 0.09(5) meV. The specific heat in zero field is expo- 
nentially activated with an activation energy A = 0.89(1) meV. A spin gap is also found through 
inelastic neutron scattering on powder samples which identify a band of magnetic excitations for 
0.8 < huj < 1.5 meV. Using sum-rules we derive an expression for the dynamic spin correlation 
function associated with non-interacting propagating triplets in a spin ladder. The van-Hove singu- 
larities of such a model are not observed in our scattering data indicating that magnetic excitations 
in Cu2(l,4-diazacycloheptane)2Cl4 are more complicated. For magnetic fields above H c i — 7.2 T 
specific heat data versus temperature show anomalies indicating a phase transition to an ordered 
state below T = 1 K. 
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I. INTRODUCTION 



As a consequences of its quantum criticality, the low-temperature properties of the Heisenberg spin-1/2 antiferro- 
magnetic chain change dramatically with small perturbations. For example, weak coupling between a macroscopic set 
of spin chains leads to Neel order at T = while modifications to the one dimensional network of spin interactions 
can induce spin-gaps and suppress Neel order. The gapped phases of one-dimensional spin-1/2 systems are simple 
examples of a qualitatively new type of magnetism where the many-body eigenstates are more conveniently described 
in terms of the singlets and triplets of near-neighbor coupled spin pairs, rather than in terms of the Zeeman states of 
individual spins. 

Theorists have established that spin gaps can be induced by bond alternation, |l],|| by frustration due to next- 
nearest neighbor interactions, |?|-^| and by coupling even numbers of chains to form spin ladders. |6[]7[| Spin-gaps have 
been observed in the quasi-one-dimensional spin- 1/2 systems Cu(N0 3 ) 2 J2.5H 2 §,§, CuGe0 3 pficj S rCu 2 Q 3 [0, 



(La,Sr,Ca)i 4 Cu240 4 i, |12j Cu 2 (l,4-Diazacycloheptane) 2 Cl4 (CuHpCl) |3|-|16[, and (VO) 2 P 2 7 Oi^fl- Apart from 
the spin-Peierls system CuGe0 3 , all of these other materials were originally thought to be spin ladders. However, 
subsequent experiments established that Cu(N0 3 ) 2 -2.5H 2 and (VO) 2 P 2 Oy actually consist of alternating spin-chains 
which are perpendicular to the putative spin-ladder. |l!^,^0| This highlights the difficulty of distinguishing between 
different Hamiltonians with spin gaps on the basis of measurements which mainly probe the magnetic density of 
states. 

This paper is an experimental study of the gapped spin-1/2 system CuHpCl. Because magnetic exchange interac- 
tions in this material are quite small it is an excellent system in which to study the field and temperature dependence 
of spin correlations in a gapped spin-1/2 system. j2lj The Cu dimer formula units in CuHpCl stack to form double 
chains. Correspondingly the magnetic interactions between Cu spins are thought to yield a spin ladder as shown in 
Fig. g with a spin-Hamiltonian of the form 

H = Ji ^2 • S 2i+ i + J 2 ^(S 2i • S 2 i +2 + S 2i _i • S 2i+ i ) + J 3 ^S 2i _i-S 2i . (1) 

i i i 

We report susceptibility, specific heat and neutron scattering experiments which determine the spin gap and magnetic 
density of states in CuHpCl. While our susceptibility and specific heat data can be interpreted in terms of a spin 
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ladder, details of the inelastic magnetic neutron scattering data are inconsistent with a simple approximation to the 
dynamic spin correlation function for a spin ladder. 



II. EXPERIMENTAL TECHNIQUES 

CuHpCl is monoclinic, with space group P1\jc and lattice constants a — 13.406(3) A, b = 11.454(2) A, c = 12.605(3) 
A and (3 = 115.01(2)°. pll The CuHpCl molecules stack to form chains in the (101) direction. The intradimer Cu-Cu 
distance is d — 3.422 A and the interdimer spacing along the chains is u = 7.00A. We have grown single crystals 
with typical masses of 5-10 mg by slow cooling of saturated methanol solutions. Crystals of this size were used for 
ac susceptibility and specific heat measurements. The zero-field ac susceptibility for T > 2 K was measured in a 
commercial SQUID magnetometer. The heat capacity was measured in a dilution refrigerator using the relaxation 
method |^2j in magnetic fields up to 9 Tesla. 

Inelastic neutron scattering measurements on a powder sample were performed using the QENS time-of-fhght 
spectrometer at the Intense Pulsed Neutron Source (IPNS) of Argonne National Laboratory. The sample consisted of 
5.75 g of powder in an annular aluminum sample can with inner diameter 0.3 cm, outer diameter 1.0 cm, and height 
12.75 cm. These dimensions were chosen to match the QENS beam size, while keeping the transmission through the 
hydrogenous sample above 50%. 

QENS is an inverse geometry time of flight spectrometer operating at a fixed final energy Ef = 3.635 meV selected 
by Bragg reflection from three horizontally-focusing pyrolytic graphite PG(002) analyzer systems. The resulting 
elastic energy resolution is SE = 0.12 meV Full Width at Half Maximum (FWHM). Each channel is also equipped 
with a cooled BeO filter which rejects neutrons with Ej > 3.7 meV through Bragg diffraction. The analyzers can 
rotate about the sample on a common rotation stage whose angular setting, v can vary from to 30°. The mean 
scattering angles, 28, of the three analyzer systems are v + 25°, 105° — v, and v + 115° respectively, and the angular 
acceptances 528 are 7.5°, 15°, and 15°. 

The magnetic neutron scattering cross section of a powder sample at wave vector transfer, Q, and energy transfer, 
Huj can be written |23| 

<|J " k '-NrA(Q)\ 2 2S(Q,u) (2) 



dridE' k Ul 2- 

where ro = 5.38 fm, N is the number of copper ions in the sample, /(Q) is the magnetic form- factor j23], and S(Q, u>) 
is the spherically averaged scattering function: 



S(Q, w)= \ d -^\ E(<W - Q a Qp)S af3 (Q, w). (3) 



The dynamic spin correlation function is given by 

lSQ/5(Q ' w) = 2^hf £ < 5 rW S r'(°) > e- lQ ' (R - R,) . (4) 

RR' 

The following corrections were applied to the neutron counts in the time-histograms associated with each of the 
three analyzer systems. First a time independent background measured for Tiw w —Ef was subtracted. Then the 
data was scaled to the shifted time-dependent count rate in a pre-sample fission monitor, and finally the data was 
converted into ftaj-histograms. This procedure yields Tiw-dependent data, I{Q, tiui), which is related to the scattering 
cross section through convolution with a resolution function as follows 

J(Q,foj)=C / dQ%cLj'K Qu (Q-Q\u-u')y-^^(Q\hLj'), (5) 

In this expression the resolution function TZq^, is assumed to be normalized to unity: 

dQ'hdu'n Qui {Q',w') = 1 (6) 



For a hydrogenous sample such as CuHpCl the elastic intensity is dominated by the incoherent nuclear scattering 
cross section which reads 
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r 

The summation here is over the nuclei associated with a single copper atom. We used the integrated elastic intensity 
to determine the spectrometer constant, NC: 

J^hdu;I(Q,tuv) 

NC = EMV (8) 

With this number we derived a normalized intensity in units suitable for comparison to theories of S(Q, to), 

The magnetic contribution to this quantity is related to the spherically averaged dynamic spin correlation function of 
Eq. (|J) as follows: 

I m (Q,hu)=2 J dQ'Mu"Jl Q UQ-Q',u-Lj')\^f(Q')\ 2 S(Q',Lj') (10) 



III. EXPERIMENTAL RESULTS 
A. Magnetic Susceptibility 

Fig. H shows the zero- field ac susceptibility of CuHpCl single crystals. The characteristic rounded peak of a ID AFM 
appears at T p = 8 K. Fitting a Curie- Weiss law to high temperature data (30 K < T < 200 K) gives a Curie- Weiss 
temperature Qcw = —6.05(1) K. Below the peak \ drops much more rapidly than for the linear S=l/2 chain with 
only NN interactions [p5|, and tends to zero at low temperatures. These data are consistent with previous powder 

] . The solid lines are a fit based on exact diagonalization of a spin-ladder model 
which will be discussed in detail below. 



susceptibility measurements [jl4fl5 



B. Specific Heat 

Figure | shows the specific heat C(H, T) of CuHpCl vs T in fixed fields up to H = 9 T. No subtraction of the lattice 
specific heat has been attempted, but as the data are strongly field-dependent and show no sign of a contribution to 
C oc T 3 out to T = 8 K, the features we observe are evidently dominated by the spin degrees of freedom. In zero 
field, C shows a broad maximum at T ss 4.5 K, below which it drops rapidly to zero. As the field is increased (H = 3 
and 5 T) the maximum in C broadens and shifts to lower T, indicating a transfer of spectral weight to lower energies. 
The qualitative low-T behavior remains the same as for H = 0, however, as may be seen from the semi- logarithmic 
plot of C vs 1/T of the H = 0, 3 and 5 T data shown in Fig. [|. All three data sets approach straight lines in the 
limit 1/T — > oo indicating that C(T) is exponentially activated. The change in slope with field indicates that the 
activation energy decreases with increasing H . This may be quantitatively related to the closure of the spin gap with 



field, as will be discussed in Sect. IV B below. 

As H is increased further, the shift of spectral weight to lower energy continues as shown in Fig. |^. As the critical 
field H c i ~ 7 T |l^Jl5[] where the gap closes is passed, the curvature of the low-T data changes, and first a shoulder 
and then a second broad peak appear at T « IK. Above H&, a cusp-like anomaly also becomes visible on top of this 
second peak. These low-T features are shown in more detail in Fig. pi where we plot C/T vs T 2 for T < 1.4 K. The 
cusp can clearly be seen at H — 8 T and H = 9 T, and below the cusp, C(T) cx T 3 . These features suggest that above 
H c \ a phase transition occurs to a field-induced ordered state. A cusp is not as clearly visible at H = 7.5 T, but below 
the pronounced knee in that data set C(T) oc T 3 , indicating that at this field the system is also ordered. In contrast, 
for H = 6.6 T, at the lowest temperatures shown C drops faster than T 3 , suggesting the existence of a small gap [ p6[ . 
The inset to Fig. || shows the H—T phase diagram of CuHpCl as determined from these data in combination with our 
earlier field-dependent susceptibility measurements H], and the magnetization data of Chaboussant et al. | fi5[ . We 
note that this phase diagram, and indeed the overall features of the susceptibility and field-dependent specific heat 
are quite similar to other gapped spin chains, such as the alternating chain Cu(NOa)2 ■ 2.5H 2 B. 
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C. Neutron Scattering 



Figure ^| shows the normalized angular average I(fiuj) of neutron scattering from CuHpCl. The normalization of 
I(fiui) is as described in Eq. ||, and the angular average was carried out over the full accessible range of scattering 
angle 21° < 28 < 153°. The intense elastic peak arises from incoherent nuclear scattering as discussed in Sect, [n] and 
the inelastic scattering for Hlo > 2 meV comes from incoherent inelastic nuclear scattering. There is however also a 
broad peak for 0.8 meV <Tllo < 1.4 meV which we now show is associated with magnetic scattering. To isolate this 
feature we measured the incoherent nuclear scattering from 5.4 g of the non-magnetic organic compound benzophenone 
(dibenzyl ketone). The material and the amount were chosen to match the total elastic cross section of CuHpCl, and 
filled approximately the same volume as the CuHpCl sample. 33 mg of boron carbide was added to the benzophenone 
to match the absorption cross-section of CuHpCl. The resulting integrated elastic intensity from the benzophenone 
sample was only 15% stronger than for CuHpCl. We applied a scale factor to the benzophenone data so that its 
integrated elastic intensity matched that of I(Q, hu) for CuHpCl. The angular-averaged data for benzophenone are 
shown as open squares in Fig. ^. As may be seen, the line shapes of the elastic peaks are indistinguishable, which 
allows us to use the benzophenone data to approximately subtract the non-magnetic contributions to the neutron 
scattering intensity from CuHpCl. A small constant term of 0.03 meV per Cu was added to ensure that the difference 
data for Tllu << —ksT averages to zero as required by detailed balance. 

The resulting difference data, I m (u>), is shown in Fig. 0(a). The only feature in these data above background is the 
peak centered at tuv = 1.15 meV. To determine whether this feature is magnetic or is due to phonons, we measured 
the scattering from the CuHpCl sample at T = 20 K and subjected those data to the same background corrections as 
the low temperature data. The result is shown in Fig. (b) . The peak broadens and a mirror image of it emerges for 
Tilo < as expected from detailed balance ]23| . This behavior establishes that the difference data shown in Fig. 0(a) 
is indeed magnetic and is a measure of the the angular averaged magnetic scattering intensity for CuHpCl: 

,■153° 

7 m (M = / <Z20/ m (Q,M, (11) 

J21° 

which may related to S(Q,u) through Eq. M. The data shows that magnetic scattering occurs only above a gap, 
A w 0.8 meV and that the bandwidth for the magnetic excitation spectrum is 0.7 meV. 

Having established the magnetic "density of states," we turn to the wave vector dependence of the magnetic scat- 
tering cross section which contains information about spatial spin correlations. In Fig. |^ we show the Q-dependence 
of the magnetic scattering integrated over the range of energies with appreciable magnetic signal: 

pl.5 meV 

L(Q) = / KdwI m {Q,hu), (12) 

J 0.8 meV 

The limited sampling density and statistics do not allow us to clearly discern any features in the data. We note only 
that as expected for magnetic scattering and not for inelastic scattering from phonons the intensity decreases with 
increasing wave vector transfer. The lines through the data in Figs. and result from a model calculation to be 
discussed below. 



IV. DISCUSSION 



A. Analysis of susceptibility data 

In their pioneering work on CuHpCl |l3f| , Chiari et al. suggested that this material may realize the Hamiltonian Tt 
of Eq. (|l|). On the basis of high temperature susceptibility data alone, they were not able to place strong constraints 
on the possible models, although they concluded that the intra-dimer and inter-dimer NN interactions J\ and J3 are 
different as indicated in Fig. |[ 

We have previously analyzed some of the data on this material in terms of the isolated dimer model ( Eq. (|l|) 
with J3 = 0, J2 = 0), the alternating near-neighbor model (J2 = 0, J3/J1 <G [0,1]), and the next-nearest-neighbor 
model without alternation (J3/J1 = 1, J2/J1 £ [0.3,0.7]). None of these provided a satisfactory fit to our powder 
susceptibility data |14| . Chaboussant et al. analyzed their powder susceptibility data in terms of both the alternating 
chain and the spin ladder (J3 = 0) models, using high temperature series expansions [Hj]. While both models fit the 
data, from structural considerations they concluded that the ladder is more likely. A comparison of their magnetization 
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data to exact diagonalization calculations for the ladder model P7J found good agreement for J2/J1 — 0.18, although 
additional improvement to the fit was found if a small ferromagnetic diagonal coupling J3 was added. 

To analyze our single crystal susceptibility data, we carried out a systematic search of the parameter space of the 
Hamiltonian in Eq. (|l|). We diagonalized chains up to length N = 14 for values of the exchange constants on a large 
grid in J2/J1 — J3/J1 space. The calculated eigenvalues for each point were used to generate x(T/Ji,H — 0) on the 
interval T/J x e [0.2,20]. Below T/Ji = 0.2, the results for N = 12 deviated from those for N = 14 by greater than 
1%, due to effects of the finite system size. By tri-linearly interpolating in this 3D mesh of (T/Ji, J 2 , J3), we could 
generate x(^) f° r arbitrary Ji, J2, and J3. The accuracy of this procedure was checked by comparing the interpolated 
x{T) between mesh points against exact diagonalization at those same values. The difference was found to be less 
than 1 part per thousand. The parameter space was constrained by the measured upper critical field H C2 = 13. 3T, 
or 18. 8K [jl5|], which occurs at g^BH c2 = max(Ji + J 3 , J-y + 2 * J 2 )) where the energy of the spinwaves from the 
ferromagnetic state change sign. 

A simultaneous fit to both the a-axis and b-axis data in Fig. || gave g a = 2.15(1), gb — 2.04(1), Ji/fcs = 13.26(4) K, 
JijkB — 2.4(3) K, and Ja/fce = 1-1(6) K. The error bars quoted correspond to the intervals over which the reduced 
X 2 deviates by less than 1/Nf ree from its minimum value. Here Nf ree is the number of data points minus the number 
of parameters in the fit. The results are shown as solid curves in the figure. The values of the g-factors in our fits 
are consistent with the values determined by EPR, Jl5| and the values for J\ and J 2 are in good agreement with the 
analysis of the magnetization data, although J3 on which our error bars are large differs in sign. 

The fit indicates that the intra rung coupling J\ exceeds the inter-rung coupling J 2 by a factor of approximately 
five. In this strong-coupling limit the zeroth order approximation to the ground state is simply a chain of non- 
interacting rung singlets and the excited states are localized triplets at energy hu) — J\. Clearly this approximation 
is too coarse since it cannot account for the finite bandwidth of the magnetic excitation spectrum which is evident in 
the neutron scattering data (Fig. ^ (a)). We therefore must take the inter-rung coupling into account, and the most 
straightforward way to do this is through a perturbation expansion. Neglecting effects of J3, the excited states then 
become propagating triplets with a dispersion relation p8[ 

Hu(q) = Ji[l + ^-cosq + 7 A 2 (3 - cos2g) 
Ji 4 Ji 




where q = Q • u is the component of wavevector transfer along the chain. For our ladder model the higher order terms 
in Eq. ( p"3| ) are very small. For example Ji(J 2 / Ji) 2 = 0.039 meV which could be compared to the energy resolution of 
our neutron scattering experiment which is SE = 0.12 meV FWHM. Therefore we approximate the dispersion relation 
with 

MQ) = Ji + J 2 cos(Q • u). (14) 
The corresponding gap in the magnetic excitation spectrum is 

A = Ji - J 2 (15) 



B. Analysis of Specific Heat Data 

To test the spin ladder model we compare its predictions to the low temperature specific heat data. If we approximate 
the dispersion relation by Eq. (|lj) and neglect interactions between propagating triplet excitations, the zero-field 
specific heat in the low temperature limit takes the form 

C=^= ( ^-) 3/2 Mexp(-A/k B T), (16) 
2V27T \k B T J V -h 

measured in J/(K mol-Cu) with n = 3 being the number of low energy modes in a Brillouin zone. The solid line 
through the zero-field data in Fig. |] is a fit to this functional form. It yields A = 0.89(1) meV, which is close to the 
value A = J\ — J 2 = 0.93(3) meV predicted by the spin ladder model. If we assume that the main effect of fields 
H < Hd on the spectrum is to Zeeman-split the triplet mode, then Eq. [lfj] should remain a good description of the 
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low-T behavior of C(H,T), although both n and A may now vary with field. The fits shown in Fig. || for H = 3 T 
and 5 T yield A = 0.55(4) meV and n = 1.3(3), and A = 0.24(4) and n = 0.6(2), respectively. This reduction of n 
with field is qualitatively consistent with spin-splitting of the triplet, as the lowest-energy band should dominate the 
specific heat in the range ksT < A(JJ)/4 where these fits were done. The gap vs field determined from the fits is 



shown in the inset of Fig. |] together with the estimated value of the gap at H = 6.6 T |26|. The observed reduction 
of A with H is consistent with the linear behavior expected for the Zeeman-splitting of the triplet. The data in the 
inset arc fit by the line shown, which has a slope of — 2.09(l)/xs, and extrapolates to A = at H = 7.2(1) T, in accord 
with H c i determined by other methods |l4|,[l5|]. 

C. Single mode approximation to S(Q,u>) 

A spin ladder Hamiltonian with J3 w thus can account for the temperature dependent susceptibility and the 
low temperature specific heat of CuHpCl. We now discuss to what extent the spin ladder model can account for our 
neutron scattering data. 

Exact diagonalization studies suggest that most of the spectral weight for magnetic excitations is associated with 
the resonant triplet mode [p9| . Therefore we use the single mode approximation: that the dynamic spin correlation 
function takes the form 

5 Q/3 (Q, tiu) = S(Q)S(Tiuj - nu(Q))6 a ,f), (17) 
Two sum-rules [po| place further constraints on 5(Q,w): 

J d 3 Q J ^^^'"(Q.w)/ J d 3 Q = S{S+l) (18) 
h 2 J udLjS^(Q,u) = < [[H,S%,Si Q ] >, (19) 

Where 

R 



From Eqs. ©, @, and (g) with J 3 = it can be shown that 



2 1 

5(Q) = ~ 



3 hu(Q) 



Ji (S 2l • Sm+i) (1 - cos(Q ■ d)) 

+J 2 «S 2j • S 2l+2 ) + (S 2l _i • S 2l+ i>) (1 - cos(Q • u)) 



(21) 



Here d and u are vectors connecting spins within a rung and along an arm of the ladder, respectively, and (Si • Sj) 
are the corresponding groundstate equal time correlation functions. Because we are interested in the strong rung 
coupling limit we neglect interdimer correlations and adopt a simplified version of Eq. (121 



S(Q) « -| M q, u) < S « • S 2«+i> (1 - cos(Q • d)) . (22) 

In this approximation the coupling along the ladder is reflected only in the one-dimensional dispersion relation Tkj(Q-u) 
of Eq. (Jlj). Requiring that S(Q,u>) defined in this way satisfies the total moment sum rule of Eq. (|lgl ) fixes 

<S 2l • S 2J+1 ) = -S(S + 1)(1 - (J 2 / Ji) 2 ) 1/2 . (23) 

With the exchange constants determined from the susceptibility measurements, this gives a model that may be 
compared to the data of Figs. and || with no adjustable parameters. 

After inserting Eq. ( p2| ) in Eq. (||) , the spherical average may be computed exactly, yielding 



S(Q h ) - _ 2 jM^Hl_^!±il \p 1 - J (dama^Q 2 - Qf) cos(Q i dcosa) 

3Q U Q~<Q f^VUl + J 2 " " (Jl - M) 

xG(Ji + J 2 -M9(fiu-(Ji- J 2 )) (24) 



6 



where Jo(x) is a Bessel function, a = 55.4° is the angle between u and d, (jl3| and the Qi — (1/u) cos^ 1 [(hu> — Ji)/ J%] 
are the solutions of Eq. (|lj) at fixed Convolution of Eq. (|24|) with the instrumental resolution and summing over 
the values of 29 used in the experiment gives the solid curve in Fig. 0(a). Several features of this curve are worth 
noting. First, the mean energy and bandwidth match the data quite well. This is to be expected since the model 
after all was chosen to account for susceptibility data whose temperature dependence controlled by these two features 
of the magnetic excitation spectrum. Second the area under the data closely matches the area under the calculated 
curve. This indicates that the inelastic peak which we have detected contains most of the magnetic neutron scattering 
to be expected from the copper ions in CuHpCl. 

The model does not, however, reproduce the shape of the data very well. The two peaks in the model correspond to 
singularities in the one-dimensional magnetic density of states at hto — J\ ± J2. The data does not show these peaks, 
being instead peaked near the center of the band. Note that inclusion of the interdimer correlations from Eq. |2l], 
which were neglected in our approximate treatment (Eq. |2^), will only make the agreement of model with data worse, 
as that produces a contribution to S{Q,uj) with a single peak at Tilj ~ J\ — Ji- We also note that the data show an 
unexpected tail of scattering intensity below the lower bound predicted by the model. This is consistent with specific 
heat data |Q that show a gap A = 0.89(1) meV, which is close to the gap seen from neutron scattering data but less 
than the gap A = J\ — J2 = 0.93(3) meV predicted by the ladder model derived from the susceptibility data. 

Fig. [sj compares the data and the model calculation of the Q-dependent energy integrated intensity (Eq. |l^). For 
CuHpCl, this is the first comparison of data which is sensitive to the geometry of the spin system to predictions of the 
ladder model. The data bear no evidence of the maximum in I(Q) predicted by the spin ladder model for Q ss 1.3A -1 , 
which raises the question of whether the ladder model is the correct description of CuHpCl. 



V. CONCLUSION 



We have reported experiments on CuHpCl which establish that this system has a gap in its magnetic excitation 
spectrum and which determine the magnetic bandwidth. Our susceptibility and specific heat data can be accounted for 
in terms of a spin-ladder model with intra-rung coupling J\ = 1.143(3) meV and two mutually frustrating inter-rung 
interactions: J2 = 0.21(3) meV and J3 = 0.09(5) meV. Not surprisingly this model reproduces the gross features of 
the magnetic density of states as probed by neutron scattering. However our neutron scattering data do not show the 
Van-Hove singularities and the Q-dependence predicted by a single mode approximation of the excitation spectrum 
in a simple spin ladder model. It is possible that a refined calculation of S{Q, uS) for a spin ladder would resolve these 
discrepancies, but it is clear that further experiments, particularly neutron scattering on deuterated single crystals, 
are required to definitively establish the spin Hamiltonian for CuHpCl. 
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FIG. 1. Schematic of chain structure of CuHpCl, indicating possible near-neighbor bond alternation (5) and 
next-near-neighbor bonds. The ellipses show how the dimer units stack. 

FIG. 2. AC susceptibility x(T) of single crystal CuHpCl with the field oriented along the a and b axes. The curves differ 
by a scale factor consistent with the g-factor anisotropy. [^| The solid lines are a fit to exact diagonalizations of the J1-J2-J3 
model. 

FIG. 3. Magnetic specific heat C(H,T) vs T for CuHpCl. At low fields C is suppressed at low temperatures because of the 
gap to magnetic excitations. Above H c i ~ 7 T, there is evidence for a phase transition to an ordered state. 

FIG. 4. Semi- logarithmic plot of C(H,T) vs 1/T in the gapped phase showing activated behavior with fits to a spin-ladder 
model as described in the text. Inset: gap A determined from fits plus an estimate of gap for H — 6.6 T pi| . 

FIG. 5. C/T vs T 2 at low temperatures, showing C oc T 3 below the temperature of the ordering anomaly at high fields. 
Inset: putative phase diagram for CuHpCl based on specific heat and magnetization data Jl4||l5| , indicating boundary between 
one- dimensional (ID) disordered and ordered (AF) phases. 

FIG. 6. Normalized, angular- averaged inelastic neutron scattering intensity for CuHpCl at T = 1.5 K (solid circles). The 
full incoherent elastic peak (small circles) is shown reduced in magnitude 1000 times. Also shown is scattering intensity for 
benzophenone (open squares), used to measure the non-magnetic background for -1.5 meV < Tioj < 2 meV. 

FIG. 7. Magnetic, angular-averaged inelastic neutron scattering intensity for CuHpCl obtained after subtracting 
non-magnetic background in Fig. 6. The reduction of the peak in (a) with increasing T identifies it as due to spin fluctuations. 
The solid line in (a) is the predicted intensity based on the spin ladder model derived from susceptibility measurements with 
no adjustable parameters. 

FIG. 8. Wavevector dependence of the energy-integrated magnetic scattering for CuHpCl at T — 1.5 K. Solid line is the 
prediction of the spin ladder model. 
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FIG. 1. Hammar, et al. 
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